Synchronization of spin-transfer oscillators driven by stimulated microwave currents 
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We have simulated the non-linear dynamics of networks of spin-transfer oscillators. The oscillators 
are magnetically uncoupled but electrically connected in series. We use a modified Landau-Lifschitz- 
Gilbert equation to describe the motion of each oscillator in the presence of the oscillations of all 
the others. We show that the oscillators of the network can be locked not only in frequency but 
also in phase. The coupling is due to the microwave components of the current induced in each 
oscillator by the oscillations in all the other oscillators. Our results show how the emitted microwave 
power of spin-transfer oscillators can be considerably enhanced by current-induced synchronization 
in an electrically connected network. We also discuss the possible application of our synchronization 
mechanism to the interpretation of the surprisingly narrow microwave spectrum in some experiments 
on a single isolated spin-transfer oscillator. 

PACS numbers: 85.75.-d,75.47.-m,75.40.Gb 



The spin transfer phenomenon, predicted by J. 
SlonczewskiQ in 1996, is now the subject of exten- 
sive experimen tal^ , E 3 , H, IM pi 0, H, @ and theoretical 
studiesJH [H 111 LI El Ea It has first been shown 
that a spin polarized current injected into a thin ferro- 
magnetic layer can switch its magnetization. This occurs 
for current densities of the order of 10 7 A. cm -2 and the 
switching can be extremely fast (< 200 ps)[l^. More re- 
cently, it has been experimentally demonstrated that, un- 
der certain conditions of applied field and current density, 
a spin polarized dc current induces a steady precession of 
the magnetization at GHz frequencies (ill ITa . Hj f . These 
steady precession effects can be obtained in Fi/NM/F 2 
standard trilayers in which a thick magnetic layer Fi 
with a fixed magnetization is used to prepare the spin 
polarized current that is injected in a free thin magnetic 
layer F 2 . The Giant Magnetoresistance effect (GMR) |2cj 
of the magnetic trilayer converts the magnetic preces- 
sion into microwave electrical signals. We will refer to 
these non-linear oscillators as "spin transfer oscillators" 
(STO). They emit at frequencies which depend on field 
and dc current, and can present very narrow frequency 
linewidths|2j]. As a consequence, they are promising 
candidates for applications in telecommunications, where 
the need for efficient, integrated and frequency agile oscil- 
lators is growing. The main drawback of the spin transfer 
oscillator is its very weak output microwave power, that 
can be optimistically estimated at - 40 dBm for a sin- 
gle oscillator. A solution to overcome this difficulty is 
to synchronize several oscillators, i.e. to force them to 
emit at a common frequency and in phase in spite of the 
intrinsic dispersion of their individual frequencies. This 
is essential for applications and, on the other hand, this 
raises complex problems which are new in spintronics and 



related to the general field of the dynamics of non-linear 
systems. 

Synchronization has been extensively studied since the 
80 's, not only because of its many potential applications 
(in physics, biology and chemistry) but also because un- 
derstanding the behavior of a large collection of non- 
linear dynamic systems is a theoretical challenge [13, El ■ 
In solid state physics, a well known example of synchro- 
nization is given by a network of Josephson junctions. An 
alternating potential takes place across a single supercon- 
ductor/insulator/superconductor junction if a dc current 
exceeding a critical current is injected through it. For an 
array of such junctions, electrically connected in series or 
in parallel, each junction emits a microwave current that 
adds to the injected dc current. When the resulting in- 
teraction exceeds a critical level, it tends to synchronize 
the oscillation of the junctions [3 El El- The theo- 
retical prediction!^ is that for N oscillators, not only 
the emitted power increases as N 2 , but the frequency 
linewidth decreases as N~ 2 . There is a definite similarity 
between networks of Josephson junctions and of STO's, 
in spite of the different equations ruling these two sys- 
tems. Recent experiments have shown that STO's can 
phase-lock (synchronize) to an external microwave cur- 
rent source [28j. Slavin et al. have analytically studied 
this case for weakly non-linear spin transfer oscillators 
|29j |. Even more recently, it has been shown experimen- 
tally that two nano-contact spin transfer oscillators can 
synchronize, but the origin of the coupling is still debated 

Bin. 

In this letter, we develop numerical simulations to 
study the synchronization of STO's. More specifically, 
for spin transfer oscillators electrically connected in sc- 
ries (or in parallel), we introduce the coupling due to 
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the microwave current induced in each oscillator by the 
oscillations of all the others. We show that, under cer- 
tain conditions for the dispersion of the frequencies, the 
GMR amplitude and the delay between the magnetic pre- 
cession and the current oscillation, synchronization can 
be obtained with an output power increasing as N 2 for a 
collection of N oscillators. 
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FIG. 1: (a) Sketch of N oscillators connected in series and 
coupled to a load Rc (throughout the paper, Rc = 50 Q). 
(b) Variation of the resistance versus the angle Qi between 
the magnetizations of Fi and F2 for an oscillator and corre- 
sponding notations. 

We first consider N oscillators of standard structure for 
spin transfer Fi(fixed)/NM/F 2 (free) connected in series 
and coupled to a dc current generator and to a resistive 
load Rc as shown on Fig^a). Our notation is displayed 
on Fijfljb). We call Kpi and Rapi the resistances of the 
oscillator i in respectively its parallel and antiparallel 
magnetic configurations. We define Roi = (RAPi+Rpi)/2, 
Ai? 4 =(R AP ,-R«)/2, (3 R = R c / (Rc+E^iRo*), and 

Pari — AR-i / (Rc+X^iR-oi)- F° r the dependence of 
the resistance Ri of the oscillator i on the angle between 
the magnetizations of Fi and F2 at time t, Qi (t) , we as- 
sume the following standard equation : 



Ri — R 



Of 



ARiCos[0i(t)] 



(1) 



The angle 0$(t) depends on the initial value of Qi at t 
= and on the variation of the current between and t. 

In first order of ^2@ARi and with the notations of 
Fig^a), a straightforward calculation leads to : 



/ = /i+5^Jij8ajm cos(9i{t)) 



(2) 



with 1\ = PrIo- Similar expressions can be found for 
oscillators connected in parallel, with different expres- 
sions for J and Pari- 

In order to study the behavior of N electrically coupled 
oscillators, we have performed simulations of the motion 
of the magnetizations nx, of the layers F2 of a collection 



considered as a macrospin without any dipolar interac- 
tion with the other m;. Its time evolution is given by a 
Landau-Lifschitz-Gilbert (LLG) equation which includes 
a standard spin transfer term proportional to the cur- 
rent. According to Eq[3 the current is the sum of the 
dc current I\ plus the coupling term ^ I\ Pari cos(9i(t)) 
and the motion equation of fhj can be written as: 



drh-j ^ — > ^ drhj 

' ~ ^lo m j x "eff + arrij x 



dt lu J WJ ' dt 

N 

+ 70 J[l + ^2 P^Ri cos(0i(t)] fhj x (fhj x M) 



(3) 



of different oscillators connected in scries. 



Each mj is 



where we have introduced the spin transfer parame- 
ter J proportional to / and expressed in field units. In 
a typical Co/Cu/Co devicefl^. a current density of 10 7 
A/cm 2 corresponds to about 10~ 2 Tesla. M is the fixed 
magnetization of all the Fi layers. The effective mag- 
netic field H e ff, is composed of an uniaxial anisotropy 
field H an , an applied magnetic field H app , and the de- 
magnetizing field Hd- All fields are in-plane (parallel to 
the direction of the fixed magnetization of Fi) except for 
the out-of-plane demagnetizing field. In the following, if 
not mentioned otherwise, we will consider the case of 10 
oscillators with H app = 0.2 T, = 1.7 T, and a Gilbert 
damping term a = 0.007 (values for Co). 

Simulations of the dynamics are performed using a 
fourth-order Runge-Kutta algorithm, with a calculation 
step of 0.5 ps. We have chosen the following random 
initial conditions: for each oscillator, the initial angles 
between the two magnetizations were randomly picked 
between and 10° for the polar angle Qi and between 
and 360° for the azimuthal angle 0j. We have checked 
that, under these conditions, the variation of the ini- 
tial conditions does not hinder synchronization. In order 
to introduce a dispersion in the behavior of the oscilla- 
tors, differences can be introduced in the anisotropy fields 
H an , demagnetizing fields or GMR ratios. We have 
checked that all these different types of dispersion give 
similar results. In this paper, we will focus on the first 
case, with the following dispersion : H an = 0.05 + (i-l)x 
0.01 in Tesla, i varying between 1 and 10. Finally, in a 
real experimental setup, the load can be connected by a 
few centimeters or meters of microwave cables to the sam- 
ple (consisting of closely connected nano-pillars). Delays 
can occur in these transmission lines between emission 
and reception of microwave currents by the oscillators 
(differences between the time t when the angles 0j(t) 
vary and the time t + r when these variations induce 
current changes). These delays are taken into account in 
the present simulations. 

We will first consider that all the oscillators have 
the same resistance R and magnetoresistance AR, 
so that we can write the coupling term in EqEl as 
JA GM R/NJ2cos{Q t (t)) with A GM r = AR/ (R -I- R c /N). 
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In this particular case, for large N, Aqmr is close to the 
value of the GMR ratio. 
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FIG. 2: Logarithm of the power versus frequency for the set of 
10 oscillators described in the text and for different coupling 
factors Aqmr- J = 0.035, and r = 5 ps. 

In FigEl we show the emitted power by the set of 10 
oscillators as a function of the frequency for different cou- 
pling parameters Aqmr- For this set, with 0.05 T < H Qn 
< 0.14 T, the dispersion of the individual frequencies is 
2.7 %, of the order of the dispersion (1.25 %) in recent ex- 
periments [3(J • The emitted power at a given frequency is 
derived by Fast Fourier Transforming the electrical power 
released in Kq- As the functions 8i(t) are known only for 
a finite number of times t, an oscillation at a well defined 
frequency appears in the Fourier transform as a peak of 
finite width. The injected dc current J is 0.035 T, and 
the delay r = 5 ps. When Aqmr, = 0, each oscillator 
oscillates at its own frequency : the frequencies are dis- 
tributed between 25.8 and 26.5 GHz approximately, and 
the total emitted power is that emitted by the sum of 
independent oscillators. For Aqmr = 0.03 and 0.05, all 
the oscillations result in a single peak. In these two cases, 
as expected, there is an increase by a factor of about 100 
in the integrated emitted power with respect to the case 
without coupling. This scaling with approximately N 2 
indicates that the N oscillators are locked not only in 
frequency but also in phase, as it will be discussed in 
more detail below. In Fig|3 we can also notice a general 
upward shift of the frequency as the coupling increases. 

In Fig|3 we consider the evolution with time of the 
trajectories of 100 oscillators. The bias conditions (J 
= 0.035 T, t = 5 ps) are similar to the previous case 
with H Qn picked randomly between 0.05 and 0.1 T for 
each oscillator (see scale on FigEJl, and Aqmr = 0.03. 
The black curve corresponds to the two symmetrical final 
trajectories. By looking at the position of the oscillators 
at different times, we see they are turning in phase (small 
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FIG. 3: Motion of 100 oscillators on their trajectory at t = 
12, 18 and 30 ns : the phase of the oscillators is locked. 



bounded phase shift), with the fastest oscillator opening 
the way. 
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FIG. 4: Frequency versus injected dc current J for oscilla- 
tor 1 and Agmr = (no coupling). The insets show the 
trajectories of the magnetization fh (M x ,M y ,M z ) in the two 
regimes referred to in the text. In the absence of current, the 
equilibrium along H e // corresponds to = 1, M y — M 2 = 
0. 

Experimentally, varying the coupling parameter 
Agmr means changing the GMR ratio in a controllable 
way, which might be difficult. Another way to increase 
the coupling is to increase J which, from Eq|3 enhances 
both the mean torque and the coupling between i and 
j. The variation of the frequency of oscillator 1 (H on = 
0.05 T) with J in the absence of coupling (Agmr = 0) 
is shown on Fig^J Similar results have been obtained in 
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simulations by other groups |17J. For J smaller than 0.016 
T, the frequency decreases as J increases in the regime of 
in-plane precessional trajectories of the magnetization. It 
increases for J larger than 0.016 T corresponding to the 
regime of out-of-plane orbits. In Fig0 we have plotted 
the difference in frequency, Si, between the 10 th oscilla- 
tor and the I s * as a function of J for different coupling 
parameters Aqmr', 5f = means synchronization of the 
two oscillators. The reference curve (no synchronization) 
obtained for Aqmr = is plotted in FigOJa). 
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FIG. 5: Difference in frequency between oscillator 1 (H a n = 
0.05 T) and oscillator 10 (H an = 0.14 T) as a function of the 
dc current. Three cases are considered, (a) Aqmr = 0, (b) 
A-Gmr = 0.03 and r = 5 ps, (c) Aqmr = 0.4 and r = 0.3 ns. 
The black arrows indicate synchronization (bo - fi =0). 

We first consider the curve of Fig|5f b) corresponding 
to Aqmr — 0.03 with a delay r of 5 ps. For low val- 
ues of J, the coupling is small, and the oscillators do 
not synchronize. The system is nevertheless disturbed 
by the injection of the microwave currents, as can be 
seen from the differences between Si for Aqmr. = 0.03 
and Aqmr = 0. Synchronization is reached (Si = 0, see 
arrows on Fig[S{b)) above J = 0.035 T (this is in the out- 
of-plane regime with, as shown on Fig[3fa), weaker dis- 
persion, and probably, easier synchronization). FigJSfc) 
corresponds to a situation with enhanced coupling (larger 
Aqmr)- In this case, synchronization extends to the in- 
plane precession regime (see arrow at J = 0.01 T). 

Another important parameter acting on the synchro- 
nization is the delay time r. We have studied the number 
of synchronized oscillators out of ten as a function of the 
delay t, with J = 0.035 T, A GM r = 0.1 and H an = 0.01 
+ (i-1) x 0.01 T. For values of t smaller than about 
0.25 ns, the number of synchronized oscillators varies pe- 
riodically with t, with a period of approximately 39 ps. 
As the frequency at synchronization is about 27.5 GHz 
which corresponds to a period T of 37 ps, there is a link 



between the period T of the oscillators and the period of 
the oscillatory variation as a function of r. The oscilla- 
tors can be synchronized more easily when they receive 
at time t microwave currents emitted at time t-nT (n = 
integer), since at t-nT they were in the same state as at 
t. [32} . This periodic behavior of synchronization ver- 
sus delay, with a period corresponding to the oscillation 
frequency, has been already predicted in other systems 
ruled by different sets of equations jU H^- The best 
conditions are for r between 0.2 and 0.75 ns and the pro- 
portion of synchronized oscillators decreases again above 
0.8 ns. 

In conclusion, we have shown that it is possible to syn- 
chronize a network of spin-transfer oscillators by simply 
connecting them electrically in series to a load (similar 
effects can be expected for oscillators in parallel). The 
synchronization depends on the dispersion of the indi- 
vidual frequencies, on the coupling parameters and the 
delay time r. Under certain conditions, the synchroniza- 
tion can be complete. In this case, the output power of 
N oscillators turns out to scale with N 2 . We have also 
shown that, for synchronized oscillators, the frequency 
as well as the emitted power are strongly dependent on 
the coupling factor Aqmr, related to the GMR ratio. 
These results are of interest for obtaining an enhanced 
microwave generation with networks of spin-transfer os- 
cillators. They also show that magnetic devices can be 
synchronized in the same way (from the coupling mech- 
anism point of view) as in the model system represented 
by a network of Josephson junctions, but with two de- 
grees of freedom (polar and azimuthal angles) instead of 
one (phase). 

We finally point out that the synchronization mech- 
anism by microwave current components we have dis- 
cussed for networks could also be important in the in- 
terpretation of the properties of a single spin-transfer os- 
cillator (pillars or point contacts). The microwave spec- 
trum of some isolated oscillators is surprisingly narrow, 
in contrast with the inhomogeneous broadening predicted 
by simulations based on micromagnetic models of ferro- 
magnetic dotsj^. However, from our results, introduc- 
ing the coupling between different parts of the dot due 
to the microwave component of the total current could 
synchronize these different parts. Such synchronization 
effects could thus explain less chaotic oscillations than 
predicted and account for the narrow linewidth of the 
microwave spectra. 
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